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Strain effects in monolayer Iron-Chalcogenide superconductors 
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The successful fabrication of one monolayer FeSe on SrTiOs represented a real breakthrough in 
searching for high-Tc Fe-based superconductors (RefP). Motivated by this important discovery, we 
studied the effects of tensile strain on one monolayer and bulk iron-chalcogenide superconductors 
(FeSe and FeTe), showing that it produces important magnetic and electronic changes in the systems. 
We found that the magnetic ground state of bulk and monolayer FeSe is the block-checkerboard 
phase, which turns into the collinear stripe phase under in plane tensile strain. FeTe, in both bulk 
and monolayer phases, shows two magnetic transitions upon increasing the tensile strain: from 
bicollinear in the ground state to block-checkerboard ending up to the collinear antiferromagnetic 
phase which could bring it in the superconducting state. Finally, the study of the mechanical 
properties of both FeSe and FeTe monolayers reveals their enormous tensile strain limits and opens 
the possibility to grow them on different substrates. 


Introduction 

Since the discovery of high critical temperature su¬ 
perconductivity in fluorine doped iron arsenide LaFeAsO 
(Tc =26 the attention of the scientific community 
has been captured by Fe-based superconductors (FeSCs). 

Among all FeSCs families known, the iron-chalcogenide 
(FeCh) (FeSe, FeTe, FeSe^^Tei-^^) family has a fundamen¬ 
tal importance for the understanding of the structural, 
electronic, magnetic and superconducting properties of 
FeSC^. In fact, thanks to their simple crystal structure, 
being composed by one layer of square Fe lattice tetra- 
hedrally bonded to Ch (see Fig0, FeCh can be thought 
as a fundamental building block of this class of materials 
and the perfect test case to study the superconducting 
mechanisms in these compounds. 

FeSe superconducts at ambient conditions with 
Tc=8 K, probably due to presence of Se vacancieJ^. On 
the contrary, FeTe is not superconductor, however, par¬ 
tial substitution of Te with Se allows superconductivity 
to emerge, reaching Tc ^15 K at the optimal substitution 
of 50%. The origin of the superconducting phase is com¬ 
monly attributed to the developing of spin fluctuations 
with wave-vector Q = (tt, tt) (associated with a collinear 
antiferromagnetic stripe phase (AFMl) in which ferro- 
magnetically ordered nearest-neighbouring Fe stripes are 
aligned perpendicularly to the antiferromagnetically or¬ 
dered ones, see Figjl|P^. In fact, despite the structural 
and electronic similarities with FeSe, FeTe shows a bi¬ 
collinear antiferromagnetic phase (AFM2]P. 

Structural properties are found to be strongly linked 
to the superconducting properties in these materials. In¬ 
deed, a strong enhancement of Tc up to 37 K is observed 
in FeSe at pressure near 9 GPsP. On the other hand, 
FeSe thin films under tensile strain of about 3.7% show 
the suppression of superconducting phas^. The role of 
the pressure and strain has been discovered to be of fun¬ 
damental importance also in FeTe. Not only pressure 
does not induce superconductivity in FeTe, but favours a 


ferromagnetic phas^. However, it was shown that FeTe 
thin films become superconductor at 13 K under tensile 
strain conditional^. 

A turning point in the exploration of FeSGs was re¬ 
cently provided by the discovery that low dimensionality 
raises dramatically the superconducting critical tempera¬ 
ture of FeS^. Recent experiments, in fact, show exciting 
properties of monolayer (ML) FeSe, with respect to the 
bulk phase. In particular, Meissner effect was observed 
at 20 K in few layer FeSe grown on SrTiOs (STO) sub¬ 
strate (lattice constant 3.90 A) El. Peng et a/.EI reported 
an enhancement of the superconducting critical temper¬ 
ature in the extremely expanded FeSe monolayer grown 
on Nb doped STO (lattice constant 3.99 A) measuring 
a superconducting gap up to 70 K. There is direct ob¬ 
servation of zero resistivity under 40 K in ex-situ exper- 
iment^i^ and ^ 109 K in the recent in-situ experiment 
with 4-point prob^^^. 

All these results represent rather solid proofs that a su¬ 
perconducting phase exists indeed in FeSe ML. 

The Fermi surface (FS) of m onolay er FeSe, measured 
using in situ ARPES techniqu^EM ^ function 

of the film thickness grown on STCP^^ reveals that it 
is characterized by only electron FS at M point of the 
Brillouin Zone (BZ), apparently losing the hole FS at F. 
Surprisingly, this is topologically different from the FS of 
the bilayer (and bulk) system, which shows the typical 
FeSG’s FS composed by hole cylinders at the BZ center 
and electron features at the BZ corner^. 

At the moment various mechanisms have been pro¬ 
posed to explain the high superconducting transition 
temperature of the FeSe ML. In particular, the electron- 
phonon coupling with oxygen optical phonons in STO 
which can induce a small momentum pairing further in¬ 
creasing the spin-fluctuation coupling^. The role of the 
substrate was also emphasized in Ref.El here the authors 
suggest that the electron-phonon coupling in FeSe is en¬ 
hanced by the presence of the substrate imposing both 
structural and magnetic constrains. 
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FIG. 1: Crystal structure of iron chalcogenides and magnetic 
phases considered. Upper panels: Side (left) and top (right) 
view of FeCh monolayer crystal. Large spheres indicate Fe 
atoms while smaller ones Ch atoms. The unit cell of the bulk 
phase is drawn as thin line and the directions of in-plane (a 
and b) and out-of-plane (c) lattice vectors are highlighted. 
Lower panels, from left to right: AFMO, AFMl, AFM2 and 
AFM3 magnetic phases (see text) with the corresponding in 
plane unit cell. 


Indeed, in these particular two dimensional (2D) con¬ 
ditions, an additional role can be played by the substrate 
and thus by epitaxial strain which may, also, suppress 
the magneto-elastic coupling. 

Considered the reported sensitivity of the physical and 
superc ondu cting properties of FeCh to structural con- 
ditionPISl, 

in this paper we studied by first principles 
density functional theory (DFT) the effect of strain on 
the structural, magnetic and electronic properties of FeSe 
and FeTe in both ML and bulk phases with the aim of 
understanding the origin of the enhancement (or of the 
suppression) of the superconducting phase. 


I. COMPUTATIONAL DETAILS 

Calculations were performed u sing the Vienna Ab- 
initio Simulation Package (VASP)P^^, using the Gen¬ 
eralized Gradient Approximation (GGA)P^ to the 
exchange-correlation energy. Van der Waals interactions, 
not negligible in FeCh system^^, were included as disper¬ 
sive term, using the DFT-D2 Grimme’s semi-empirical 


correctiorP^^^. 

We studied the strain effects considering different mag¬ 
netic phases (see below). In order to simulate the struc¬ 
tural effects of the STO substrate (which constrains the 
in-plane lattice constants), we neglected the in-plane 
magneto-elastic distortion. Once the ground state ge¬ 
ometry for each system was determined including mag¬ 
netic distortions with inequivalent a and b directions, 
the systems were constrained to in-plane lattice constants 
a = (|a| + |b|)/2: this is assumed to be the reference sys¬ 
tem in the un — strained condition. Strain effects were 
studied calculating the total energy as a function of the 
applied in-plane strain (varying the a in-plane lattice con¬ 
stant) considering different competing magnetic phases 
(schematically shown in FigQ: the checkerboard order 
of the iron spin (called AFMO), the collinear stripe phase 
(AFMl), the bi-collinear one (AFM2) and the recently 
proposed block-checkerboard order (AFM3)P^. The mag¬ 
netic order along the out-of-plane axis is considered fer¬ 
romagnetic for all the bulk phases. We used Projected 
Augmented-Wave pseudopotential^ for all the atomic 
species involved, with an energy cutoff of 350 eV. Inte¬ 
gration over BZ was performed using different uniform 
Monkhorst and Pack grid^ depending on the lattice: 
16 X 16 X 10 for the tetragonal (a x a x c) non mag¬ 
netic (NM) phase (containing two Fe and two Ch atoms), 
8 X 8 X 10 for magnetic cell (2a x 2a x c) commensu¬ 
rate to the AFMO, AFMl and AFM2 magnetic phases. 
The AFM3 phase was simulated using y/2a x 2v^a x c 
unit cell using a Monkhorst and Pack grid of 8 x 4 x 10. 
For the ML systems we considered a supercell including 
^ 10 A of vacuum. Using these parameters the total en¬ 
ergy of the considered systems are converged to less than 
0.1 meV/atom. 

The calculations of band structures and FS’s have been 
performed using the NM electronic state^^. 


II. RESULTS AND DISCUSSION 

The total energy of FeSe ML as a function of the in¬ 
plane lattice constant for the lowest energy magnetic 
structures are shown in Fig. [^a). At variance with what 
was expected in line with the collinear stripe ground state 
of pnictides, the ground state magnetic order is the block- 
checkerboard with equilibrium structural parameters re¬ 
ported in Tab. |Tj Similar with what found in the AFMl 
phase, the magneto-elastic coupling causes the inequiva¬ 
lence of a and b directions. 

The magnetic moment of Fe atoms in the monolayer 
phase is 0.2 {Ib higher with respect to the bulk system 
(see Tableand changes slightly (less than 0.1 iib) de¬ 
pending on the magnetic order. 

FigJ^a) shows that the AFMl phase is the more sta¬ 
ble phase under applied in-plane strain, which induces 
a magnetic phase transition at a 3.9 A. Consider¬ 
ing that the superconducting pairing could be mediated 
by (tt, tt) spin-fluctuations, free-standing stoichiometric 
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FIG. 2: Total energy per atom as a function of the in-plane 
lattice parameter for ML FeSe (panel (a)) and FeTe (panel 
(b)). The theoretical points are htted with a cubic function. 
The dotted vertical line indicates the NbiSTO lattice param¬ 
eter (see text). 


AFMO o 
AFMl • 
AFM3 A 


•X 5)^ 

\ \ 


NbiSrTiOs, 


p A. 

0 '' 


.•. 

. A .A-. 


-©. 


.5 3.6 3.7 3.8 3.9 4.0 4.1 

a (A) 



'(b)' 

AFMl • 



AFM2 ■ 



AFM3 A NbiSrTiOg 



A-',.''' 



/ •''' 


\ ^ 

/■ F 


V\ 

> X' 


\ "x 



..^ 



- 


point and two (quas i ove rlapping) electron-like FS’s at 
the BZ corner (Fig. 3(a)). Stretching the FeSe ML at 
the NbiSTO lattice parameter (Fig. 3(b)[ ) fills completely 
the dx 2 _y 2 hole-type band near F and makes the corre¬ 
sponding FS to vanish. At the same time, the electron 
FS’s change their shape and shrink as a result of electron 
counting. ARPES measurements, however, seem to find 
a different scenariJ^ only electron sheets at the BZ cor¬ 
ner are found with no evidences of hole-type FS at the 
BZ center. 

Although the calculated electronic properties do not 
completely agree with available experiments, we indeed 
predict that strain induces a magnetic phase transition 
and a subsequent topological transition in the FS of the 
ML phase of FeSe. 



FIG. 3: Fermi surface of FeSe at equilibrium (a) and strained 
at NbiSTO lattice parameter (b). 


FeSe monolayer with a (tt, 7t/2) magnetic order is not ex¬ 
pected to be superconductor, however it may support a 
superconducting phase if strained at the STO lattice con¬ 
stant in agreement with experimental evidenceJI^. The 
AFMl phase increases its stability under applied strain 
and becomes the ground state at the NbiSTO lattice con¬ 
stant (see FigJ^a)) where experiments report the highest 
critical temperatur^I^. 



FeSe 

FeTe 

|a| (A) 

3.75 

4.04 

|b| (A) 

3.72 

3.51 

f^ch (A) 

1.44 

1.79 

\m\ (hb) 

2.1 

2.5 


TABLE 1: Equilibrium structural parameters of FeSe and 
FeTe ML in their respective ground state (|a| and |b| lat¬ 
tice constants and the height of the chalcogenide, hch) and Fe 
magnetic moment, \m\. 

The variation of about 6% of the in-plane lattice con¬ 
stant (NbiSTO) determines a strong effect on the FS, as 
reported in Eig. This is mainly due to the contraction 
of the Se height from the Ee plane, hch (FigQ, which 
moves from 1.44 A (equilibrium geometry) to 1.34 A at 
the in-plane lattice constant of NbiSTO. In the free¬ 
standing ML, there are three hole-type ES’s at the F 


We now consider the monolayer phase of FeTe, a sys¬ 
tem not yet experimentally realized. As already dis¬ 
cussed, bulk FeTe is not superconducting, however the 
effect of the in-plane strain has been shown to induce 
a superconducting phase with a critical temperature of 
13 In principle, a monolayer of FeTe can result 

more interesting than the FeSe monolayer; in fact, if su¬ 
perconductivity is mediated by spin-fluctuations, FeTe is 
expected to have a stronger pairin^^. 

The ground state of FeTe ML is found to be AFM2 
with lattice parameters reported in Tab. |lj While the 
magnetic induced distortion of the square lattice is mini¬ 
mal in FeSe, in FeTe it is sensible. In addition, the mag¬ 
netic moment of Fe atoms is larger in the ground state 
(AFM2) of FeTe than in FeSe (see Table and shows a 
reduction of 0.3/j.b in the AFMl phase. From Fig. |^b) 
we note that two magnetic phase transitions take place 
as a function of the tensile straini from AFM2 to AFM3 
at a > 3.80 A and from AFM3 to AFMl at a > 4.05 A. 

The FS’s shape of FeTe ML changes from thn equi¬ 
librium (Fig, 4(a)) to the strained phase (Fig. |4(b) ). In 
this last geometry, the FS shows two hole pockets at 
F and two (quasi overlapping) electron features at M 
nearly nested with Q = (77,77). Therefore the magnetic 
transition towards the AFMl phase with this ordering 
vector allows the pairing between hole and electron FS’s 
sheets which was hindered in the bicollinear and block- 
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checkerboard phases. Based on these structural, mag¬ 
netic and electronic properties we propose strained FeTe 
ML as a good candidate for interesting superconducting 
properties, potentially better than FeSe ML. 



FIG. 4: Fermi surface for ] 
and strained at 4.17 A in- 



(b) a = 4.17 A. 

ML at equilibrium (panel a) 
: lattice parameter (panel b). 


The results obtained in the monolayer systems help 
the interpretation of the effects of in-plane strain on the 
superconducting properties of FeCh bulk systems. In¬ 
deed, for FeSe and FeTe thin films a different behaviour 
under tensile strain was experimentally reported. FeSe 
grown on STO shows suppression of the superconduct¬ 
ing phase observed on unstrained compound^, while, 
bulk FeTe becomes superconductor upon the same tensile 
strain conditional^. 

Total energy calculations performed constraining the 
in-plane lattice parameter and fully relaxing all the oth¬ 
ers {i.e. out-of-plane axis and internal positions) show 
that the magnetic ground state of bulk FeSe at equi¬ 
librium is the AFM3 order, and it becomes AFMl at 
the NbiSTO in-plane lattice constant. Thus, supp ression 
of the superconducting state in strained FeSe^^^ can be 
interpreted considering that strain causes a topological 
transition leading to a complete closure of the elect ron F S 
cylinders at M (present at the equilibrium, see Fig |5(a)[ ). 
At a > 3 .99 A the FS becomes three dimensional (see 
Fig. |5(b) ), strongly reducing the pairing. 

On the contrary, we find that FeTe bulk has a bi- 
collinear (AFM2) magnetic ground state in agreement 
with experiments^. However, a transition towards the 
AFM3 phase takes place at a 3.8 A; upon further in¬ 
crease of the tensile strain above a > 4.10 A, the sys¬ 
tem develops an AFMl ground state. In this last phase, 
bulk FeTe preserves both hole and electron sheets al¬ 
though with sensible differences with respect to the equi¬ 
librium phase (FigJ^. On the basis of these results, 
the experimental evidences reported by Han et al^ find 
now a natural interpretation: strain causes AFMl spin- 
fluctuations which couple the hole and electron FS sheets 
via the AFMl Q vector. 

In order to possibly guide new experiments towards 
the use of alternative substrates to properly tune strain 
conditions in both bulk and monolayer systems, we in¬ 
vestigated the mechanical properties of FeSe and FeTe 



(a) a = 3.73 A. 



(b) a = 3.99 A. 

FIG. 5: Fermi surface for FeSe bulk at equilibrium (panel (a)) 
and strained at NbiSTO lattice parameter (panel (b)). The 
F-point is the center of the first BZ shown as thin lines. 


ML to obtain information on the stiffness and breaking 
strength of these compounds. Starting from the equilib¬ 
rium structure, we studied the stress-strain relation for 
FeSe and FeTe ML and report the results in Fig. [7l Here, 
the tensile strain is defined as 5 = [Sxx + ^yy )with 
^xx = {a' — a)/a and Cyy = {b' — 6)/P^, being a and b 
the lattice parameters along x and y directions respec- 
tivel3@. 

The stress-strain curve for FeSe and FeTe shows a max¬ 
imum at the critical strain of 5 = 0.33 and 5 = 0.39, 
respectively. Although the corresponding maximal stress 
= 11-9 GPa and = 9.1 GPa) is very low 

compared to other monolayers as graphene or dichalco- 
genides {i.e. M 0 S 2 ), these results demonstrate that both 
compounds can support relevant straim comparable with 
what recently found for phosphorene^. This can open 
the way to future growth of superconducting chalco- 
genides on different substrates with quite large mismatch. 
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III. CONCLUSIONS 



(a) a = 3.75 A. 



(b) a = 4.10 A. 


Using first-principles DFT simulations we calculated 
the mechanical, electronic and magnetic properties of 
FeCh in 2D and bulk geometries under in-plane strain. 
Our predictions show that: 

{i) The magnetic ground state of monolayer and bulk 
FeSe is the block-checkerboard. Application of tensile 
strain produces a magnetic transition towards the AFMl 
phase with the closure of only one FS at the F point in 
FeSe ML. 

In the bulk geometry, strain induces a topological tran¬ 
sition of the FS at M, that can be at the origin of the 
suppression of the superconducting state. 

{ii) If FeTe ML could be strained above 9%, it will 
achieve all the requirements needed to observe a su¬ 
perconducting transition: AFMl magnetic ground state 
with hole and electron FS’s sheets nested with the proper 
ordering wave vector. The same effect is predicted also 
in the bulk phase as well. 

{in) Both FeSe and FeTe ML have large enough me¬ 
chanical flexibility to support critical strains up to 5 > 
30%. 


FIG. 6: Fermi surface for FeTe bulk at equilibrium (panel (a)) 
and strained at 4.10 A lattice parameter (panel (b)). The 
F-point is the center of the first BZ shown as thin lines. 



FIG. 7: Stress-strain relation for FeSe and FeTe ML. 
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